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Abstract. We report the detection of the extremely red object (ERO), HR 10 (7 - JsT ~ 6.5, 2=1.44), at 4.9 
and 6.1 /im (rest- frame) with ISOCAM, the mid-infrared (MIR) camera onboard the Infrared Space Observatory 
(ISO). HR 10 is the first ERO spectroscopically identified to be associated with an ultra-luminous IR galaxy 
(ULIG) detected in the radio, MIR and sub-millimeter. The rest-frame spectral energy distribution (SED) of HR 
10 is amazingly similar to the one of Arp 220, scaled by a factor 3.8±1.3. The corresponding 8-1000 fim luminosity 
(~ 7x 10^^ h'^Q Z/q) translates into a star formation rate of about 1200 h^Q Mq yr~^ if HR 10 is mostly powered 
by star formation. We address the key issue of the origin of the powerful luminosity of HR 10, i.e. starburst versus 
active galactic nucleus (AGN), by using the similarity with its closeby clone, Arp 220. 
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1. Introduction 

HR 10 (or ERO J164502+4626.4, Dey et al. 1999), is the 
first and presently only Extremely Red Object (ERO, usu- 
ally defined as galaxies with I — K > 4) known to be as- 
sociated with the class of ultra-luminous infrared galaxies 
(ULIGs, LiR= L(8-1000/im)> lO^^ Lq). It was detected 
by Hu & Ridgway (1994, hereafter HR) together with an- 
other ERO (HR 14 or ERO J164457-I-4626.0) in the field 
of the QSO PC 1643+4631A (z= 3.79). HR initially sug- 
gested that both galaxies with extreme colors {I — K > 6) 
could be distant ellipticals lying at z ~ 2 — 3. More gen- 
erally, deep near IR (NIR) surveys indicate that EROs 
present the same clustering properties (Daddi et al. 2000, 
McCarthy et al. 2001) and surface brightness distribution 
(Moriondo, Cimatti & Daddi 2000) as elliptical galaxies. 
But dusty starbursts being potential progenitors of local 
ellipticals, they may also show similar clustering proper- 
ties and some local examples, as NGC 7252 (Hibbard et al. 
1994) or Arp 220 (Scoville et al. 2000), already show a de 
Vaucouleurs luminosity profile typical of ellipticals. High 
resolution NIR imagery and spectroscopy with the Keck 
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telescopes (Graham & Dey 1996) revealed that HR 10 was 
a moderately distant (z = 1.44) galaxy with an asym- 
metric morphology and Ha in emission. Another evidence 
against HR 10 being an early-type galaxy is its strong 
sub-millimeter luminosity as measured with SCUBA at 
the JCMT (Cimatti et al. 1998, Dey et al. 1999). The de- 
tection in HR 10 of a large CO luminosity, hence molecular 
hydrogen mass, was recently reported by Andreani et al. 
(2000) and presented as evidence favoring a star forma- 
tion origin for the bulk of the IR luminosity rather than 
an AGN. However, a large gas mass may not only feed 
star formation but also gas accretion onto an AGN (see 
Papadopoulos et al. 2001). Another test for the presence 
of an AGN in a dusty galaxy is the warm over cold dust, 
i.e mid IR (MIR, 3-40 /im) over far IR (FIR, 40-300 ^im), 
luminosity ratio as well as the shape of the MIR spectrum. 
We present the first detection of HR 10 in two MIR bands 
of ISOCAM onboard ISO corresponding to the rest-frame 
3.3-6.1 and 4.9-7.4 /xm wavelength ranges. A comparison 
to local galaxies SEDs is presented. We show that the SED 
of HR 10 is a scaled version of that of the closest ULIG, 
Arp 220. We discuss the origin of the luminosity of HR 10 
on the basis of this similarity. 
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Fig. 1. MIR contours over /-band image (WHT archives); N: up, E: left, a) 15 /im contours at 2.5 to 5-cr, with a 0.5-ct 
step, b) I2fj.m contours at S/N=1.6, 1.9, 2.2, 2.8, 3.5, 4.1, 6, 9, 11. 



Throughout this paper, we will assume Ho= 70 km s ^ 
Mpc"-^, rimattor = 0.3 and r^A = 0.7. For this cosmology, 
the luminosity distance of HR 10 {z = 1.44) is 10373 Mpc. 

2. Observations 

Deep images of a field of 13' x 13' centered on the QSO PC 
1643+4631A (z = 3.79) were obtained with the ISOCAM 
broadband filters LW3 (12-18/im, centered at 15 /im) and 
LWIO (8-15 /xm, centered at 12 /im) for a total integration 
time of 16 and 13 minutes respectively. The 15 /im image 
results from the coaddition of two overlapping mosaics 
slightly rotated one with respect to the other resulting in 
an improved spatial resolution (2" pixels). The 12 /im im- 
age being made of one single mosaic of 6" pixels suffers 
from a lower spatial resolution. A given position of the sky 
was observed by 18 (20) different pixels at 15 /im (12 /im), 
resulting in a better flat-fielding and correction of cosmic 
ray impacts. The data reduction and source extraction 
was done with PRETI (Pattern REcognition Technique 
for ISOCAM data, Starck et al. 1999). The gaussian noise 
(photon plus detector noise) in the 15 fim (12 /im) im- 
age is 58 /iJy (31/iJy). This field was also observed with 
ISOPHOT onboard ISO at 90 and 170 /im (P.I. R.Ivison). 
The total exposure time per sky position is 11.5 and 6.6 
minutes respectively. HR 10 is not detected with a 5-(T 
upper limit of ~ 200 mJy in both bands. 

On Fig. |l|, the MIR contours are overlayed on the I- 
band image of a field containing both HR 10 and PC 
1643+4631A. Both galaxies are detected in each MIR 
band with the flux densities given in Table |l|. The un- 
certainty on the flux density related to the correction for 
the transient behavior of the detectors and to the pixel 
size relative to the point spread function is given in a 
second line in italics. While the 15 /im contours are cen- 
tered on the optical positions, the 12 /im contours of HR 
10 present an offset of 4". The astrometry was first cal- 
culated over the whole 13' field using several stars from 
the US Naval Observatory Catalog which were detected in 



the MIR. A relative astrometric correction was then ap- 
plyied to check the position of the MIR detection at the 
optical position of HR 10 by using a set of six objects in- 
cluded in a field of 2'x2' centered on HR 10 and including 
PC 1643-I-4631A. The probability of a chance association 
with an optical object with an /-band magnitude lower 
than / within a distance d was estimated using the follow- 
ing formula assuming a Poissonian distribution of sources: 
P = 1 - exp[-n(/)7rd2]. The offset of 4" for the 12 /im de- 
tection is attributed a probability of a chance association 
with the optical position of HR 10 of 20 %. This is not sur- 
prising because of the lower spatial resolution of the 12 /im 
image (6" pixels). However the probability of a random 
association of the 12 /im source within 4" of the 15 /im de- 
tection of HR 10 is equal to 0.1 % (n[0.2 mJy at 15 /i,m] = 
0.8±0.1 ' ^). In order to quantify the risk of detecting 
"ghost" sources produced by cosmic ray residuals in the 
MIR images, we performed Monte-Carlo simulations by 
inserting fake sources in real datacubes. Down to a flux 
density limit of S^, ~ 100 /iJy, we find that about 6 and 
12% of the sources are false in the 15 and 12 /im images 
respectively. However, we were able to reduce the fraction 
of ghost sources to zero in the 15/im image by requiring 
that a source be detected in both mosaics (before coaddi- 
tion) at a lower detection level. Only above 0.3 mJy do we 
reach the zero probability of a false detection in the 12 /im 
image, but again the probability of a random association 
of a spurious 12 /im source with the 15 /im detection of 
HR 10 is negligible. The upper limits on HR 14 [I — K ~ 
6.2, see Table |l|) imply that it cannot be a ULIG unless 
it is more distant than z ~ 1.8. PC 1643-I-4631A is one of 
the most distant sources detected in the MIR [z = 3.79 
corresponds to an age of the universe of only 1.6 Gyr) and 
belongs to the class of ULIGs, with rest-frame luminosi- 
ties of: vL^[i.l ± 0.63 /im] = (2.3 ± 0.5) x lO^^ h^^ Lq 
and i/i,,[2.5±0.84/im] = (1.2 ±0.8) x lO^^ /j-2 
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Fig. 2. Rest-frame SED of HR 10 (filled circles, Dey et al. 1999) including our MIR data, a) comparison to the SED 
of Arp 220 normalized by a factor 3.8. Upper limits from IRAS (Xohs— 25, 60fim) and ISOPHOT (90, 170 jim). Plain 
line: combined spectra from ISO-LWS (43-185 iim, Fischer et al. 1999) and ISOCAM-CVF (5-18 ^m, Charmandaris et 
al. 1999). Dashed line: interpolation between observed magnitudes: UV (RC3), BI (Surace et al. 2000), RcJHKL (Kim 
1995), IRAS (FSC), ISOPHOT (Klaas et al. 1997), sub-millimeter (Rigopoulou et al. 1996), radio (Anantharamaiah 
et al. 2000). b) comparison to the mean SEDs of a "high reddening starburst" (dashed line, filled circles) and a Seyfert 
2 (dash-dotted line, empty circles) from Schmitt et al. (1997), and to Mrk 231 (dotted line, Ivison et al. 1998). 
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Table 1. MIR flux densities at 15 and 12 /im, signal-to- 
noise ratios and probability of chance association. 



3. Discussion: Nature of HR 10 

The least-square fit of the SED of HR 10 using the 
SED of Arp 220 (Fig. |^a) illustrates the amazing sim- 
ilarity of both spectra. The 8-1000 /im luminosity of 
Arp 220 is Ljr = 1.8 x 10^^ hj^ Lq, for a dis- 
tance of d= 81.8 h^fl Mpc, corrected from the Virgo in- 
fall {cz— 5439-1-206 km s^^). The resulting normaliza- 
tion factor of 3.8 ± 1.3 (l-cr) implies an IR luminosity of 
LiR = (6.8±2.3) X 1012 fQj. jjj, -^Q j^^Yie dust was 

mostly heated by massive young stars, then this IR lumi- 
nosity would translate into a star formation rate {SFR) 
of about (1170 ± 396) /if^f Mq yr-\ for a 10-100 Myr 
continuous burst, solar abundance and a Salpeter IMF 
{SFR [M0yr-i] = 1.72 x 10"!° Ljr, Kennicutt 1998). 

A consistent normalization factor also applies for the 
molecular gas mass of HR 10, estimated from its CO lu- 
minosity (Mh2= 1.2x10" Mq, Andreani et al. 2000), 
which is 2.9 times larger than in Arp 220 (AfH2[Arp 
220]= 4.1 xlO^" Mq, ScoviUe, Yun & Bryant 1997). If the 
molecular gas mass over dust mass ratio of HR 10 is the 
same as the one measured in the center of Arp 220, i.e. 



Mgas/Mdust ~ 400 (ScoviUe, Yun & Bryant 1997), then we 
obtain Mjust [HRIO] ~ 4 x 10^ Hjq Mq , in agreement with 
the one estimated by Dey et al. (1999, ~ 7 x 10* hj^ Mq) 
using a black body fit to the FIR part of the SED with 
Tdust - 40 K. 

The SED of Arp 220 is limited to wavelengths above 
the U-band (3550 A) which corresponds to 8670 A in the 
observed frame of HR 10. In Fig. we have fitted the 
SED of HR 10 with the model STARDUST2 (Chanial et 
al., in preparation) from the observed B-band (1800 A in 
the rest-frame) to the radio. If we assume that there is 
only one component responsible for both the optical and 
IR light in HR 10, then the best-fit is obtained for a vi- 
sual extinction of Ay = 3 mag (± 0.2 mag, extinction 
law A{X)/Ay from Calzetti et al. 2000), with an age of 
0.2 Gyr at the time of the observation and a character- 
istic timescale for the starburst of t— 95 Myr {SFR= 
Mgas{t)/T). The optical- UV part of the spectrum alone 
can be fitted with an Ay— 2-2.8, an age of 2-4 Gyr and 
T= 0.4-0.9 Gyr, resuhing in Ljr - 8 x 10^° /ifg^ Lq which 
does not fit the MIR and sub-millimeter part of the SED. 
However, detailed studies of local luminous IR galaxies 
have shown that the region from which the bulk of the IR 
luminosity arises can contribute very weakly to the opti- 
cal light, e.g. the Antennae galaxies (Mirabel et al. 1998). 
In such a scenario, a second component with a larger Ay 
would contribute dominantly to the IR regime and nearly 
not to the optical part. In both cases, we obtain an IR lu- 
minosity and SFR within the error bars associated with 
the fit using Arp 220. The main difference between the 
model and Arp 220 's SED is that the FIR over MIR lumi- 
nosity ratio of Arp 220 is higher. The fit of Arp 220's SED 
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requires an Ay > 30 mag for the component responsible 
for the bulk of the IR light. 

The key issue remains to determine whether the bulk 
of the luminosity of HR 10 is powered by star forma- 
tion or accretion by a black hole. Two other EROs have 
been spectroscopically identified and detected in the MIR 
and radio (ISO J1324-2016, Pierre et al. 2001, and ERO 
J164023+4644, Smith et al. 2001). The upper limits estab- 
lished for their bolometric IR luminosities in the absence 
of a detection in the FIR or sub-millimeter are consistent 
with their belonging to the class of ULIGs too. NIR spec- 
troscopy in both galaxies favor the presence of an AGN. 
In the case of HR 10, the width of the observed H^ -I- [Nil] 
emission feature or the [Nil] /Ha ratio are too uncertain 
to be used as an indication of the presence of an AGN 
(Graham & Dey 1996). The Ha line is strongly affected 
by dust extinction (S'Fi?[Ha]~ 80 Mq yr"^ only). 

Most authors favor the starburst hypothesis as a dom- 
inant source of energy in Arp 220: its MIR spectrum up 
to 40 fxia shows no evidence for high ionization lines ex- 
pected for AGNs (Sturm et al. 1996), its "IR excess" 
(LiR,/L[Lya] ^ 24, Anantharamaiah et al. 2000) is much 
lower than for AGNs (~ 45 — 65) and typical of starburst 
galaxies (~ 12 — 45, Genzel et al. 1998), its radio emis- 
sion is produced by several compact sources (Smith et al. 
1998), the ratio of aromatic features over MIR continuum 
(Genzel et al. 1998) and the slope of the MIR continuum 
(Laurent et al. 2000) are typical of starbursts. More re- 
cently, Haas et al. (2001) suggested that the MIR lumi- 
nosity of Arp 220 could be underestimated because of dust 
extinction in the MIR and that after dereddening, its FIR 
over MIR luminosity ratio would be closer to the one for 
AGNs. But the dereddening factor varies by a factor five 
depending on the dust geometry assumed. Finally, the flat 
2-10 keV hard X-ray spectrum of Arp 220 implies that in 
order to be mostly powered by an AGN, it would need to 
be Compton thick with a column density larger than 10^*^ 
cm~^ (Iwasawa et al. 2001). 

In Fig. ||b, the SED of HR 10 is compared to the mean 
SEDs of a Seyfert 2 (Sy2) and of a starburst with high red- 
dening (SBH), from Schmitt et al. (1997), normalized to 
the 15 /im luminosity density of HR 10. The 1-a error bars 
corresponding to the distribution of the 15 SBH and 15 
Sy2 galaxies are shown in the radio and at Arost= 100 /im. 
Sub-millimeter data arc missing in these SEDs but they 
should decrease in luminosity density above 100 /im as in 
local galaxies, e.g. the Seyfert 1 Mrk 231 {d— 180.9 h^Q 
Mpc, LiR = 3.5 X 10^^ Lq, dotted line). Without 
MIR data, Mrk 231 was a candidate template for HR 10 
(see Dey et al. 1999). Although not as red as HR 10, the 
SBH fits the MIR-FIR region of HR 10 while the FIR 
luminosity of the Sy2 is about ten times fainter. 

Finally, even if the presence of a combination of an 
AGN and a starburst is still an option for both HR 10 
and Arp 220, most studies favor a dominant contribution 
from star formation to their IR luminosities, implying that 
HR 10 presents the largest SFR known at present. This 
study clearly shows the need for direct MIR and FIR/sub- 



millimiter observations of distant dusty galaxies to im- 
prove our understanding of a population of objects which 
plays a major role in galaxy evolution (Chary & Elbaz 
2001), hence emphasizes the importance of the next gener- 
ation IR satellites to come, i.e. SIRTF, FIRST and NGST. 
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